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We report on an extensive experimental study of intersubband excitations in the THz range 
arising from the couphng between a quantum well and a zero-dimensional metal-metal microcavities. 
Because of the conceptual simplicity of the resonators we obtain an extremely predictable and 
controllable system to investigate light-matter interaction. The experimental data is modelled by 
combining a quantum mechanical approach with an effective medium electromagnetic simulation 
that allows us to take into account the losses of the system. By comparing our modelling with the 
data we are able to retrieve microscopic information, such as the electronic populations on different 
subbands as a function of the temperature. Our modelling approach sets the base of a designer tool 
for intersubband light-matter coupled systems. 

PACS numbers: 



I. INTRODUCTION 

The ability of engineering semiconductor devices at the 
nanometer scale, by modifying their quantum proper- 
ties, has been an essential ingredient that has allowed 
the continuous impressive technological development of 
electronics and optoelectronics. A prominent example 
is provided by the emergence of new laser sources in the 
Mid-IR and THz regions, the so called Quantum Cascade 
Lasers In these structures the electronic transport, 
under the application of an electrical bias is carefully de- 
signed in order to obtain a redistribution of the carriers 
on the different energy levels and finally population in- 
version between subbands of multiple-coupled semicon- 
ductor quantum wells. 

The engineering of the light-matter interaction in an 
intersubband (ISB) system resides somehow on the op- 
posite limit, as it can be enhanced by achieving a large 
electronic population on the fundamental subband of a 
quantum well inserted in a microcavity with a very small 
volume. In this case the electromagnetic response of the 
system is dominated by collective plasmonic excitations 
of the 2D electron gaz [2] . When combined with an opti- 
cal microcavity, these systems enter the strong coupling 
regime, in which the interaction between the electromag- 
netic and the electronic collective modes yield new mixed 
intersubband polariton states [s^. The splitting between 
the mixed states, also called Rabi splitting, is a direct 
measure of the light-matter coupling strength. The latter 
can be increased through doping by increasing the num- 
ber of electrons participating to the interaction. This 
enabled to reach the ultra-strong coupling regime, where 
the Rabi splitting becomes comparable with the energy 
of the ISB plasmon mode This regime has been in- 



tensively studied and reported by several groups, even at 
THz frequencies and up to room temperature [5|-[8|. 

The aim of this work is to interpret the spectroscopic 
features of the polariton states as a function of the tem- 
perature in order to extract microscopic information on 
intersubband transitions. In this respect we will show 
that measurements of the Rabi splitting, 2Qr^ permit 
to follow the electronic distribution on the different sub- 
bands. Indeed, 2Qii has a straightforward dependence 
on the subband population and the corresponding tran- 
sition oscillator strengths. This is clearly illustrated in 
the formula of the Rabi splitting for a single ISB excita- 
tion between the fundamental and second subband, cou- 
pled with the lowest order TMq mode of a double-metal 
waveguide (ol. fioj: 
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Here e is the electron charge and is the vacuum 
permittivity. The effective electron mass m* and the 
dielectric constant of the semiconductor core e are well 
known material parameters. The parameter Lcav is the 
cavity thickness, A/'qw is the number of quantum wells 
contained in the cavity, N1/N2 is the electronic popula- 
tion of the fundamental/second subband, and /12 is the 
transition oscillator strength that depends on the overlap 
between wave- functions. Eq.([T]) indicates that the pop- 
ulations, Ni and N2 can be inferred from the study of 
the splitting 2Qji as a function of the temperature, by 
applying the Fermi-Dirac statistics. 

Our study is conducted with OD double-metal disper- 
sionless square-patch resonators, as described in FigTO 
which operate on standing wave-like TMq modes [ll|. 
The conceptual simplicity of this geometry reduces con- 
siderably the photonic degrees of freedom, allowing the 
direct study of the electronic component of the polari- 
ton states. Indeed, these cavities feature a flat photonic 
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FIG. 1: (a) Schematic of an individual metal- metal patch cav- 
ity with thickness Lcav, containing 25 GaAs/AlGaAs quan- 
tum wells. The structure is probed in reflectivity measure- 
ments, with an incident angle 6, in both TM (p-) and TE(o-) 
polarization, (b) Microscope picture of the top surface of the 
sample, containing an array of patch cavities. The size of the 
patches is s, and the period of the array is p. 



dispersion, and the resonant frequency is independent 
from the incident angle of the probing beam. Moreover, 
the electric field that couples to the intersubband po- 
larization is homogeneous along the growth axis of the 
quantum wells. The only parameter that governs the 
coupling is the cavity mode detuning from the ISB reso- 
nance, which in turn is deterministically set by the size 
s of the patch 

Although a more complex fabrication procedure, the 
metalHc patch cavities have several advantages with re- 
spect to the multi-pass absorption geometry that is com- 
monly used to study the intersubband transitions [12|. 
First of all they allow performing intersubband transi- 
tion spectroscopy with light propagating normal to the 
wafer surface and therefore with a much larger surface 
exposed to the impinging. Moreover, in this geometry 
the light does not traverse the semiconductor substrate 
and avoids all the spurious absorptions in the bulk ma- 
terial. We should nevertheless carry in mind that Eq.([T]) 
describes an ideal system, which does not take into ac- 
count all the details of a typical experimental sample. 
To correctly deduce the electronic distribution from the 
experimentally measured Rabi splitting one needs to con- 
sider several other complications, which are typically not 
included in an purely energy preserving Hamiltonian de- 
scription: 

i) An important aspect is the dissipation in the system, 
which causes finite linewidths both for the cavity modes 
and intersubband resonances. Therefore the measured 
Rabi splitting is actually reduced with respect to values 
predicted by Eq.© [13]. 

ii) In highly doped samples, several ISB resonances can 
be active at the same time. Their interaction causes a 
redistribution of the oscillator strength [l^. Therefore 
we need a multi-transition generalization of Eq.([T]). 



iii) For the THz region, where large quantum wells are 
employed, the Hartree potential induces important cor- 
rection to the heterostructure potential. Therefore en- 
ergy levels and wave functions should be obtained by self- 
consistent Schrodinger-Poisson calculation. This should 
be taken into account for the correct application of the 
Fermi-Dirac statistics. 

iv) The presence of metal-dielectric junctions gives rise 
to a charge transfer, which changes the carrier concen- 
tration of the quantum wells close to the junction [l5|. 
This effect is important for the correct estimate of the 
number of charged quantum wells A/'qw- 

The structures that we investigate, already reported in 
Ref. [l6|, are detailed in the subsequent sections. Our 
study is based on experimental results that were not pre- 
sented in that reference, and on an analysis that takes 
into account the afore mentioned points i)-iv). 

The paper is organized as follows. In part [IT] we 
present the experimental system and we briefly review 
the photonic confinement in our structures. Part IIIII 
deals with an extensive study of the intersubband po- 
laritons at different temperatures and for different cavity 
detuning with respect to the ISB resonances. We show 
how the information collected in these experiments, in 
combination with a self-consistent simulation of the het- 
erostructure potential allows to infer the electronic pop- 
ulation of different subbands as a function of the temper- 
ature. In particular, we interpret our absorption data at 
THz frequency as a function of the temperature [iG^ , and 
show that they actually display the transfer of oscillator 
strength Q. 



II. EXPERIMENTAL SYSTEM 
A. Microcavity 

The photonic structure employed in our experiments is 
summarized in Fig(T]that sketches a single square-patch 
double- metal microcavity containing 25 GaAs/AlGaAs 
modulation doped quantum wells. The characterization 
of the quantum wells will be presented in details in the 
next section. The sample contains a dense array of such 
micro-cavities, as illustrated in Fig(T]^b), in which the 
period of the array p is kept close to the square size s, 
so to optimize the coupling efficiency with the incident 
radiation [11]. The overall thickness of the core region 
between the two metals is Lcav =1.5/im. 

The structure is tested by varying the incident angles 
(as defined in Figdja)) and the electric field polarization 
of the incident beam. Characteristic reflectivity spectra 
for a cavity made of an undoped GaAs thick layer (dotted 
line) and for a cavity with quantum wells (solid line) are 
shown in Figl2fa). Data are taken at room temperature 
(300 K) and for an angle of incidence of ^ = 45°. For both 
samples the patch sizes are equal to s = 12 jam. These 
spectra display dips at the energies of the resonant cavity 
modes of the structure. For both structures we observe 
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FIG. 2: (a) Typical room temperature reflectivity spectra 
for two samples with an array of patches with s —12fim. 
Dashed line: undoped reference sample. Solid line: sample 
with doped quantum wells, but with the same overall thick- 
ness (l^cav — 

1.5/im) as the reference, (b) Resonance frequency 
of the doped sample, measured at 300 K, as a function of the 
inverse square size multiplied by the speed of light (c/s). K 
denotes the order of the resonant mode, and the grey area 
corresponds to the optical phonon band. 



two resonances which correspond to the excitation of the 
first two TMo-like standing wave modes [lH, around 3 
THz and 6 THz. In the case of the cavity with quantum 
wells, at room temperature, there is no clear signature of 
the ISB resonances. However the effect the electron gas 
is still visible in the data of Figl2fa) through the con- 
siderable enlargement of the cavity peaks for the doped 
sample. Indeed, at high temperature, the electrons pro- 
vided by the Si-donors can scatter in the confined and 
continuum states of the heterostructure by giving rise to 
sizeable free carrier losses. 

In Figl2fb) we summarize the resonant frequencies of 
measured doped structure at T = 300 K as a function of 
the inverse square size multiplied by the speed of light: 
c/s. Measurements for incident angles = 45° (circles) 
and = 10° (triangles) are both reported. At low fre- 
quencies we recover the typical linear dispersion of the 
TMo guided mode, vk = Kc/ueffS^ with K the order 
of the lateral standing wave and Uef / the effective index. 
At higher frequencies, the dispersion becomes non-linear 
due to coupling with the optical phonon of the semicon- 
ductor, which results in a strong increase of the effective 
index. The dispersion is independent from the incident 
angle, however close to normal incidence {0 = 10°) only 
the resonances with odd ^are excited due to the selec- 
tion rule of the grating [11]. The position of the reso- 
nances is well described by a modal method model for a 
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FIG. 3: Typical conduction band profile of the periodic quan- 
tum well stack, including self-consistent Hartree corrections. 
We have highlighted a single quantum well together with the 
moduli squared of the first 4 wavefunctions, as well as the 
Fermi- level energy, Ep- 



line ar g rating |17j , based on the scattering matrix formal- 
ism jl8|. Indeed, the modes of the square patches of the 
structure, that couple better to the incident beam, have 
a dipolar charge distribution [ll|, and these modes are 
identical to the resonances of the double-metal lamellar 
grating structure. We will therefore use this formalism 
for the data modelling that is described in part IIIII 



B. Quantum well slab 

The semiconductor region is a GaAs/Alo.isGao.rsAs 
heterostructure made of 32nm thick quantum wells 
(QWs) and 20nm thick barriers. X-ray analyses indi- 
cate a slightly lower average Al content in the barriers 
12% than the nominal 15%. The structure is J-doped in 
the each barrier, 5nm away from the adjacent quantum 
well, with a nominal value of 2 x lO^^cm"^ of the sheet 
carrier concentration. 

In Figl3] we show a typical conduction band profile 
obtained with Poisson-Schrodinger solver with periodic 
boundary conditions, at low temperature (T =7 K), in- 
cluding the Hartree corrections to the quantum well po- 
tential. The envelope wave-functions of the confined sub- 
bands of the quantum well are also indicated. Details on 
these simulations are given in Appendix [A] and in part 
lllli The absorption of the quantum well is measured in 
a multi-pass configuration with 45° polished facet and a 
gold layer deposited on top, as indicated in Fig H] (inset). 
The results of the experiment are summarized in FiglH 
where two absorption peaks are clearly visible from 7 K 
to 150 K. These peaks correspond to the 1^2 and 2^3 
transitions of the quantum well of Figl3l The contour 
plot of Figm^a) is obtained by merging multi-pass trans- 
mission spectra performed at different temperatures. In 
the contour plot the red color scale corresponds to a low 
transmission signal and therefore strong ISB absorption. 
In the lower part of the figure we report the spectrum at 
T =7 K, where the 3.5 THz peak corresponds to the 1^2 
transition, and the 4.8 THz peak is the 2^3 transition. 

As expected from Fermi-Dirac statistics, the intensity 
of the two peaks have an opposite dependence on temper- 
ature: the first 1^2 transition is activated at low tem- 
perature as electrons accumulate on the ground state of 
the quantum well and the population difference A^i — N2 
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FIG. 4: (a) Multipass absorption of the quantum well 
medium, measured as a function of the temperature and fre- 
quency. The two features around 3 THz and 5 THz corre- 
spond respectively to the 1^2 and 2^3 transitions, (b) The 
absorption spectrum at T =7 K. The inset illustrates the ex- 
perimental configuration. 

is increased. Moreover, a strong blue shift is visible for 
the 1^2 transition, due to the depolarization effect [2]. 
The transition 2^3 follows an opposite behaviour, since 
the subband 2 is progressively populated as the temper- 
ature is risen. The fact that the 2^3 transition is visible 
even at low temperature indicates that the Fermi level 
lies above the second subband edge, meaning that the 
doping is higher than the nominal value. Further, we 
shall see that the ISB absorption of the 2^3 transition is 
increased due to the oscillator strength transfer between 
the first and the second ISB excitation of the quantum 
well [14]. Namely, this phenomenon will be elucidated 
and quantified from the study of the splitting induced by 
the light-matter interaction on each individual transition. 



C. Polariton dot measurements 

When the QW medium is inserted in a OD microcavity 
the strong interaction between the ISB excitation and 
the cavity mode gives rise to two polariton states. In our 
system the spectroscopic features are a function of the 
size of the resonator 5, which sets the detuning between 
the cavity mode and the ISB excitations. 

The polariton splitting due to the coupling between 
matter excitations and the first order mode of our 0- 
dimensional cavities is very insensitive to the experimen- 
tal conditions. This can be seen in the series of different 
experimental conditions reported in the data of Figl5]and 
Figini These experiments are performed at low temper- 
ature, T =7 K, where the 1^2 ISB excitation is coupled 
with the fundamental K = 1 cavity mode. This mode is 
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FIG. 5: (a) Experimental configuration for reflectivity mea- 
surements with variable polarization of the incident beam, 
(b) Low temperature (T =7 K) reflectivity spectra performed 
at an incident angle = 45°. Solid curve: TE linear po- 
larization. Dashed curve: unpolarized light. The size of the 
patches is s=ll.l/im. 
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FIG. 6: (a) Experimental configuration for reflectivity mea- 
surements with different incident angles, (b) Low temper- 
ature (T =7 K) reflectivity spectra performed with un- 
polarized light, for an angle = 10° (solid line) and = 45° 
(dashed line). The size of the patches is s=12.5//m. 



twofold degenerate into TMioo and TMqio modes, each 
oscillating along the two sides of the square stripe (along 
the X- and y- direction respectively). Since the two modes 
have identical coupling with the ISB transition, the po- 
lariton splitting is independent from the particular polar- 
ization state of the incident beam, as described in FigEl 
where we compare two spectra obtained with unpolarised 
(dotted hue) and polarized (solid line) light. If the inci- 
dent beam is ^/-polarized, then only the TMqio mode is 
excited. The resulting polariton reflectance spectrum, 
obtained for an angle of incidence ^ = 45° is shown in 
solid line in FiglS] The size of the patch in this exper- 
iment is 5=ll.l/im. The two polariton peaks are indi- 
cated as = 1' and K = V . The same experiment re- 
peated with unpolarized light (dashed line) gives rise to 
practically an identical spectrum (dashed curve in FigEJ. 
As a matter of fact for = 45° both TMioo and TMqio 
are excited and therefore the polariton states K = V and 
K = 1" can be similarly excited for the polarization in 
the X- and y- directions. Contrarily, one can see that the 
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FIG. 7: Frequency - temperature contour plots of the reflectivity for cavities with different sizes s. The resonant frequencies of 
the cavities (indicated by arrows) have been chosen so that the microcavity modes are tuned between the 1^2 and 2^3 ISB 
transitions. The dashed lines indicate the values of the ISB resonances measured at high temperature, T > 140 K, according 
to the data in FiglH 



second order TM200 mode (indicated as = 2) has a 
different dependence and appears only if a component of 
the incident electric field in the x- direction is present. 
This is due to a symmetry selection rule of the even res- 
onances [11]. 

In Figl6]we report the measurements performed with 
unpolarized light, but at different incident angles: = 
10° (solid curve) and = 45° (dashed curve). The patch 
size is s=12.5/im in this case. Once more we can see that 
the polaritons {K = V and K = V^) are unaffected by the 
angle of incidence, since they belong to the coupling with 
the first order mode. The K = 2 resonance disappears at 
normal incidence due to its symmetry selection rule [ll[ . 

The above properties of the polariton states, that at- 
test their OD nature, stem directly from the properties of 
the K = 1 resonance, which is dispersionless (indepen- 
dent from the angle of the incidence) and insensitive to 
the polarization because of the twofold degeneracy aris- 
ing from the symmetry of the structure. These properties 
are imprinted on the ISB polarization through the strong 
interaction with light. We therefore refer to these struc- 
tures as "polariton dots" |16i] . 

Having reduced the light degrees of freedom as much as 
possible, we can now concentrate on the study of the po- 
lariton characteristics as a function of their matter part, 
i.e. of the ISB excitation. To this end we have measured 
the temperature evolution of the polaritons for different 
cavity resonances, as illustrated in Figd The four pan- 
els in Fig|7] represent the reflectivity contour plots as a 
function of the temperature for cavities, in which the first 
order {K = 1) resonance is varied between the 1-^2 and 
2^3 ISB transitions. As an indication, the dashed lines 
in FiglTlreport the values of the ISB resonances measured 
at high temperature, T > 140 K (for graphical simplicity 



we do not consider here the temperature dependence of 
the ISB resonances). Note that these values are closest 
to the ones that one would expect for the same QW with- 
out charge effects. In the color code, red corresponds to 
low reflectance, and blue corresponds to close to unity 
reflection. The cavity resonance is indicated with an ar- 
row. The first panel (cavity with s=12.5/im) is resonant 
with the 1^2 transition, while the fourth (5=7.8/im) is 
resonant with the 2^3 transition. The first panel is per- 
formed at 45° angle of incidence, while the others are 
performed at 10°. This is why in the first panel it is 
visible also the second order resonance {K = 2) at 6 
THz. One can see that as the temperature is lowered 
below 150 K electrons accumulate preferentially on the 
first subband, making the ratio N1/N2 >> 1. The 1^2 
transition is thus activated resulting in a progressively 
increasing of the polariton splitting for the 5=12.5/im 
cavity or, a growing absorption dip for the off-resonant 
cases. The 2^3 transition is most active in the tempera- 
ture range (50 K < T < 150 K), as qualitatively expected 
from the Fermi-Dirac statistics. 

The main advantage of the microcavity measurements 
with respect to multi-pass measurements is the type of 
spectroscopic data that is related to the population dif- 
ferences between the QW subbands. In a multi-pass ex- 
periment, such as the one depicted in Fig HI if a linear 
absorption regime is assumed, the population difference 
Ni — N2 is related to the integral of the absorption peak 
in the spectra [12]. The peak area determination requires 
a careful baseline calibration, which can become difficult 
in the THz region owe to diffraction effects due to the 
large wavelengths, as well as an increased free carrier 
and phonon absorption. In a microcavity, the popula- 
tion difference A^i — N2 can be inferred precisely from the 
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polariton splitting (the diflFerence between the two polari- 
ton frequencies) observed in the reflectivity experiment. 
Moreover, both methods require a careful determination 
of the linewidth of the ISB resonances. 

The impact of the dissipation effects on the spectra, 
such as free carrier absorption and ISB linewidth enlarge- 
ment is visible in the data of Fig|71 For high temperatures 
the cavity resonances appear to have rather low quality 
factors, whereas the quality factors improve at low tem- 
perature, becoming comparable with the values measured 
for an undoped structure. This is clearly visible, e.g. in 
the case of the 5=7.8/im cavity, where much narrower 
cavity resonance is recovered after the 2^3 transition is 
depleted at very low temperature: the FWHM is 1.2 THz 
at 200 K and 0.79 THz at 7 K. This behavior confirms 
the comparison between the doped and undoped refer- 
ence (FiglSfa)). We attribute it mainly to the enlarge- 
ment of the ISB resonances, that couple to the cavity 
mode, which is also visible in the multi-pass absorption 
spectra (FigH]). Moreover, at high temperature the free 
carrier absorption in the heterostructure also increases, 
as part of the electronic population can be thermally ex- 
cited in the high levels of the well and the 3D continuum 
of the heterostructure. 



To carry out this approach, we implement self- 
consistent Poisson-Schrodinger coupled equations to cal- 
culate the potential changes associated with the field 
produced by the static charges. As a result, we obtain 
the envelope wave- functions (piiz) and energy levels hwi 
produced by the electronic confinement of the quantum 
wells, as wells as the Hartree correction. If the total sheet 
concentration Ntot is known, we can then obtain the pop- 
ulation of each subband A^^ through Fermi-Dirac statistic. 
We can then use this information in order to formulate 
the light-matter coupling Hamiltonian [10]. The inter- 
action of the transition i ^ j with light can be fully 
described introducing the function ^ij(^), which repre- 
sents the oscillating electronic current between levels i 
and j: 



dz 



(l)j{z) 



dz 



(2) 



In the following, to simplify the notation, we use a 
single greek index a to indicate the transition i ^ j. 
For instance, the frequencies of the ISB transitions are 
uja = ojj — uJi . The current functions of Eq.(j2j) contain 
all the relevant quantities of the problem. Namely, the 
oscillator strength of the transition a is: 



III. DATA EXPLOITATION 

To ascertain the above qualitative discussion we have 
modelled the experimental data for cavities with different 
dimensions at different temperatures. The data mod- 
elling is based on an effective medium approach, com- 
bined with numerical simulations of the optical response 
of the double metal structures. This approach allows us, 
on one hand, to consider the situation of several occupied 
subbands, and on the other, to take into account the fi- 
nite linewidths of both the microcavity mode and the ISB 
resonances. The ohmic losses in the metallic layers are 
taken into account by considering a complex Drude con- 
stant for the metal [HI, [HI The effective medium constant 
can be obtained from the quantum Hamiltonian describ- 
ing the light-matter interaction [lOj. In the following we 
start by recalling the quantum approach. 



A. Plasma Hamiltonian and Effective medium 
approach 

The description of light-matter interaction using a 
Hamiltonian approach allows to take into account all the 
microscopic details of the electronic confinement. More- 
over, the effective dielectric function of the quantum well 
media can be obtained from the characteristic polyno- 
mial of the quantum Hamiltonian. The quantum model 
allows therefore to obtain a dielectric function that is 
both directly applicable for electromagnetic simulations 
and encapsulates the microscopic details of the system. 



and the effective thickness Lry of the transition is: 



e^{z)dz (4) 



Further, the effective length Eq.(|4]) allows us to define 
the plasma frequency of the transition uopa'- 



(5) 



where AN^ = Ni — Nj is the population difference 
(surface density) between levels i and j. With the above 
definitions Eq.(|5]) coincides with the plasma frequency 
introduced by Ando, Folwer and Stern Finally, we 
introduce the frequency of the intersubband plasmon uj^ 
which takes into account the depolarization effect: 



^Pa 



(6) 



These quantities permit to write the bosonized Hamil- 
tonian describing the interaction of the intersubband 
plasmons with the cavity mode, as well as the coupling 
between the plasmons associated with different transi- 
tions [lOj : 
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^ = ^ huJaPaPa + huJca^{a^ a + 1/2) 

+z ^ hfla{a^ - a){pl + Pa) 

a 

+ ^^^c./3(p1 +Pa)(p^ +P/3) (7) 

Here Pq, and pl^ are bosonic operators describing the in- 
tersubband plasmons, a and are creation/annihilation 
operators of the cavity mode, and cjcav is the mode fre- 
quency. The strength of the hght-matter couphng is pro- 
vided by the coefficient Q^' 

= ^.f^f^ (8) 

where = La/Lcav is the overlap between the mi- 
crocavity and the intersubband plasmon. The strength 
of the inter-plasmon coupling is provided by the coeffi- 
cient Eai3' 

The Hamiltonian ([7j) describes the coupling between 
set of quantum oscillators, one of which is the electro- 
magnetic mode. The plasmonic oscillators are also cou- 
pled between each other through Coulomb dipole-dipole 
interactions, which are expressed in the last term of the 
Hamiltonian ([7|). As described in Appendix [Bit he matter 
part of the Hamiltonian ([7]) can be diagonalised numeri- 
cally thus obtaining a set of independent plasmon modes 
each separately coupled with the light field: 

H = Y^ hWjP\Pj + nc^cav(a^a + 1/2) 
J 

+^E^^/^(«'-«)(^i+^^) (10) 

Here the Pj and Pj are the new plasmonic operators 
with frequencies Wj. The coefficients Rj play the role of 
effective plasma frequencies and characterize completely 
the coupling of the J*^ mode to the electromagnetic field. 
As an illustration, let us consider the case of a single in- 
tersubband transition, J = a. Then Rj = R^ becomes: 

Rl=^Pa^fU^ (11) 

This equation shows that the coefficient R^ contains 
altogether the information of the population differences 
AA/'q, between the subband levels, through the plasma 



frequency Eq. (j5j) , as well as the overlap between the cav- 
ity modes and the intersubband polarization. 

When several intersubband transitions are present, 
they interact through the last term of the Hamiltonian 
([7]). After the Hamiltonian is diagonalized into its form 
([To]), this interaction is also contained in the coefficients 
Rj. Moreover, owe to the plasmon-plasmon coupling, 
the coefficients Rj are different from the ones that one 
would obtain for a single uncoupled plasmon (Eq. (pT]) ). 
This is precisely the phenomenon of oscillator strength 
transfer, as explained further. 

As shown in details in Appendix [Bj the characteristic 
polynomial of the Hamiltonian (p^Oj) allows to define the 
effective medium dielectric constant of the system: 

^ = 1 + V ?^ fl2) 

where we have introduced the phenomenological 
linewidths, Fj, of the plasmon modes. This expression is 
similar to the plasmon pole approximation for 3D plas- 
mons [20]. Note that, when rotating wave approximation 
is assumed, Eq.([T2j) leads to a Lorenzian lineshape with 
full width at half-maximum equal to Vj. Once more, 
Eq. (p!2|) underlines the importance of the effective plasma 
frequencies Rj. Indeed, from the semi-classical point 
of view, R^j is the weight of the plasmon pole Wj of 
the inverse dielectric function, quantifying its oscillator 
strength. With this respect, the quantity R^j can be con- 
sidered as a generalization of the concept of oscillator 
strength in the case of ISB plasmons. 

Note that, in the case of a multiple quantum well sys- 
tem as considered here, the Hamiltonian ([7]) contains the 
couplings between plasmons from spatialy different quan- 
tum wells, as well as between plasmons from the same 
quantum well. However, in the case where we can neglect 
the overlap between wavefunctions from different quan- 
tum wells, the plasmon coupling coefficients ([9]) vanish, 
and the remaining non-zero couplings are only those be- 
tween plasmons from the same well. If the A^qw quantum 
wells are identical, we can consider a Hamiltonian of the 
form of Eq. (p!Q|) . where the operators Pj describe the 
normal plasmon modes of a single quantum well, but the 
plasma oscillator strengths are replaced by R^j x Nqyj. 
This fact is briefly discussed in Appendix [Bl 

To take into account the intrinsic anisotropy of the 
intersubband system and the in-plane free carrier ab- 
sorption, the semiconductor slab is modeled by a diag- 
onal dielectric tensor with components £z = SQf^{uo) and 
= = ^ll(^) [21]. The in-plane dielectric function 
e\\{uo) describes the Drude absorption arising from the 
free movement of the carriers in the plane (x,?/)[22]. We 
have also considered the contribution of the LO phonons 
in the dielectric function, as in Ref.0. 

The knowledge of the effective medium dielectric con- 
stant expressed in Eq. ([T2]) gives the electromagnetic re- 
sponse of the ISB system, and therefore it is a central 
quantity for the interpretation of our data. In particular 
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FIG. 8: Block diagram of the fitting procedure used to exploit 
the experimental data. Top: fitting the experimental spectra 
(reflectivity and absorption) with the electromagnetic model, 
based on the dielectric function of Ea. (|12|) provides a phe- 
nomenological estimate for the quantities ^ and Vj 

for each temperature. Bottom: the quantities Wj are 
computed from the microscopic quantum model, using the 
total number of electrons Ntot in the well. For each temper- 
ature T we retain the value of Ntot which best reproduce the 
experimental values R^f^ and Wj^^. 
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FIG. 9: Conduction band profile and charge density simula- 
tions of the system, taking into account the metal semicon- 
ductor interfaces f23'l. 



wavefunctions (j)i{z). 

We should bear in mind several underlying approxi- 
mations. First, the quantum Hamiltonians of Eq.([7j) or 
Eq. (p!Q|) describe a perfect electronic system, as they do 
not include any dissipative effects that may decrease or 
destroy the coherence of the electronic wave- functions, 
such as disorder or scattering from phonons or interface 
roughness. We approximate these effects by the inclu- 
sion of the phenomenological homogeneous broadening 
parameters Fj in Eq. (p!2]) . Moreover, the plasma Hamil- 
tonians of Eq.([7]) and Eq. (p!Q|) are valid in the limit of 
high electronic densities, where the exchange-correlation 
effect can be neglected. 



Eq. (p!2]) involves the quantities Wj, Fj and Rj describ- 
ing the collective plasmon modes. As it is often the case 
we do not know all the exact details of the system and 
therefore we regard the quantities Wj, Fj and Rj as free 
parameters that can be determined by fitting the experi- 
mental spectra. Indeed, the positions of the intersubband 
peaks Wj and their linewidths Fj can be directly de- 
rived from the multi-pass absorption experiment (FigH]). 
Moreover, the quantities Rj can be deduced from the po- 
lariton splitting in measurement such as those presented 
in FigH 

The values of Wj and Rj can also be calculated 
from the microscopic model described above, based on 
the Poisson-Schrodinger modelHng of the heterostructure 
and the diagonalization of the Hamiltonian of Eq. ([7]). 
In this approach, the critical input parameter is the num- 
ber of ionized impurities that gives rise to the electronic 
charge (sheet) density Ntot- For each temperature T, 
the sheet density Ntot is adjusted so that the calculated 
values of Wj, and Rj coincide with the values obtained 
from the fit of the experimental spectra, Wj^^^ and R^^ . 
The scheme for relating the experimental spectra to the 
microscopic computation is summarized in FiglSl Fol- 
lowing this approach we can connect the results from the 
spectroscopic measurements with the microscopic infor- 
mation on the heterostructure as a function of the tem- 
perature T. In particular we can deduce the population 
of the different subbands Ni , as well as the bare intersub- 
band energies hwi and the corresponding single particle 



B. Modelling of the experimental data 

In our structure the crucial quantity for the evalua- 
tion of the light matter interaction is the total number of 
electron participating to process. In a metal-dielectric- 
metal cavity the Fermi level pins at the metal- Ga As in- 
terface [15] and therefore some of the quantum wells in 
the immediate vicinity of the metallic layers are depleted. 
In order to estimate the number of quantum wells that 
are still charged, Nqw, we have used a freeware Poisson 
band-structure solver [235. The conduction band profile 
simulations are presented in Figl9l The simulation shows 
that, below 150 K, 2 QWs are depleted at each interface. 
Both ends of the structure are therefore modelled as un- 
doped spacers. The electronic distribution in the charged 
QWs is practically uniform, as can be seen from FiglQl 
Indeed, we model this region with an effective dielectric 
constant of the form expressed in Eq.([T2]), which assumes 
an identical average doping for all the remaining 21 doped 
QWs. To take into account the spacer regions, accord- 
ing to Eq. (pT]) the values for Rj have been multiplied by 
^21/25 = 0.92. 

The total population Ntot{T) and its redistribution on 
the different subbands as a function of the temperature 
Ni{T) can be extracted by injecting it into the Hamil- 
tonian (fTQ|) . calculating Wj and Rj (for J = 1 ^ 2 
and J = 2 ^ 3) and finally comparing these values with 
those obtained by the fit with the experimental spectra. 
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Examples of the experimental fit are provided in FiglTOl 
In Fig lTOT a) the transmission experiment without cav- 
ity (dotted curve) is compared with a transfer matrix 
simulations based on the effective dielectric constant ([12]) 
in order to extract the values of Wj and Tj. These val- 
ues are then inserted in the simulations of the reflectivity 
spectra presented in Figures [TOT b) and [TOT c), where we 
consider the particular case of cavities in resonance with 
the 1^2 and 2^3 ISB excitations respectively. From 
these simulations one can extract the effective plasma 
frequencies Ri^2 and i?2^3- Note that the reflectivity 
spectra from the cavities have intrinsically better quality 
than the multipass transmission data, especially at a high 
temperature. Indeed, in the multi-pass absorption light 
pass through a thick substrate (350/im) to probe a very 
thin active region (1.5/im). The substrate transparency is 
strongly affected by the increased optical phonon absorp- 
tion as the temperature is raised. This is clearly visible 
in the data of Fig lTOT a). where the transmission signal 
drops considerably around 7 THz. On the contrary, in 
the case of the cavity there is no substrate and the light 
field is crosses exclusively the QW region. This configu- 
ration is therefore more appropriate for testing ISB elec- 
tronic absorption and avoiding spurious effect from the 
substrate. 

The parameters Wj, and Rj obtained from this proce- 
dure have been plotted as a function of the temperature 
in Figima) and FigHU^b) (fuh circles). In Figinib) the 
full widths at half maximum of the absorption peaks, Fj, 
are also indicated as errors bars. The continuous lines are 
the results of our simulation using a Poisson-Schrodinger 
solver described in Appendix [A] together with the diag- 
onalization of the plasmon Hamiltonian ([7]). As we have 
already mentioned, we have left, for each temperature T, 
the total sheet electronic density per well as a free fitting 
parameter. In Fig(T2fa) we present Ntot as a function 
of the temperature as well as the population differences, 
{Ni — Nj){T) between the occupied subbands. For the 
fit we have used the value of 31.8nm for the quantum 
well thickness (instead if the nominal value 32nm), which 
gives us a good fit over the whole temperature range. Ac- 
cording to these results the number of electrons Ntot{T) 
captured by the quantum wells is constant below 60 K, 
and then progressively decreases as the temperature is 
raised. This behaviour can be explained by assuming 
that some of the electrons are recaptured in the delta 
doped donor region. 

In Fig JTTT b) we also report the frequencies 0Ji^2 and 
(dashed lines) calculated as the difference between 
the single particle energy states. Due to the Hartree po- 
tential in a modulation doped structure uJi^2 decreases 
as the electronic population is increased. Therefore, as 
the temperature is lowered between 7 K and 60 K, the 
energy spacing hwi^2 is reduced and the subband 2 be- 
comes populated. Note also that the collective plasma 
frequency Wi^2, that fits very well the experimental 
data, has exactly the opposite behavior with respect to 
cji^2- In the case of the 2^3 transition the depo- 



larization effect is much smaller and the collective mode 
frequency I^2^3 is close to the single particle frequency 
(^2^3. The observed slight depolarization shift for this 
transition is even lower than that predicted by the quan- 
tum Hamiltonian ([7]) that describes a perfectly coherent 
electronic system {Tj = 0). The small discrepancy arises 
from disorder effects [24]. As a matter of fact it can be 
noticed that ^2^3 — ^2^3 < r2^3, which implies that 
the correlations introduced by the depolarization effects 
are not enough to overcome the disorder, and the system 
stays closer to its single particle energy. In this case the 
effects of the inhomogeneous broadening become impor- 
tant, and models like those exposed in Ref. [25] are more 
pertinent to describe the system. 

The microscopic model used to fit the experimental 
data permits also to evidence the interaction between 
ISB plasmons from different transitions which leads to 
transfer of the oscillation strength between the first and 
excited ISB transitions [14]. Let us recall how this phe- 
nomenon is contained in the Hamiltonians ([7]) and (p^Oj) . 
The Hamiltonian ([7]) describes a set of collective modes 
at frequencies cD^. Their plasma frequencies uop^ (Eq.(j5j)) 
not only quantify their coupling strength with the cav- 
ity mode, but also describe their interaction through the 
dipole-dipole Coulomb coupling, Eq.®. By diagonaliz- 
ing the matter part of the Hamiltonian ([7]), we obtain 
the Hamiltonian (p!Q|) which contains the new indepen- 
dent normal collective modes at frequencies Wj. There- 
fore the light-matter interaction occurs between plas- 
mons originating from individual ISB plasmons coupled 
by Coulomb interaction. A side effect of the coupling 
is that, when two ISB excitations are present, as in our 
system, the first excitation transfers part of its oscillator 
strength to the second one. 

For our system this effect is not visible in the excita- 
tion frequencies, as the normal mode frequencies Wi^2 
and 1^2 ^3 remain close to the original collective mode 
frequencies cDi^2 and 002^^ (not shown). However, the 
transfer of oscillator becomes apparent by considering the 
plasma energies. This is illustrated in Fig(T2fb), where 
we compare Ri^2 and i?2^3 (continuous line) with the 
values Ri^2 ^2^3 that would be obtained for two 
completely uncoupled plasmons (dashed line). The later 
can be obtained analytically from Eq. (pT]) . We observe 
that always i?2^3 > ^2^3 and Ri^2 < Ri^2' Notice 
that at low temperature i^2^3 is almost twice ^^2^3. This 
means that the effective oscillator strength of the plas- 
mon resonance is doubled when we consider the Coulomb 
dipole-dipole interaction. This explains why 2^3 ab- 
sorption peak in the spectra of FigJIOTa) is quite high in 
spite of the relatively low population on the subband 2. 



C. Photonic dispersion and finite linewidth eff*ects 

To proof the consistency of the values of Wj and Rj 
extracted from the fits in FigITO]we have simulated the 
spectra for all the measured cavities, at four different 
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FIG. 10: Examples of electromagnetic simulations (solid curves) of the experimental spectra (dotted lines) for different tem- 
peratures, (a) Multipass transmission experiment, (b), (c) Reflectivity measurements at ^ = 45° for arrays of patches with 
dimensions s =12.5/im and s =7.8/im respectively. 
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for the two plasmon resonances, as extracted from experi- 
ments (dots) or computed from the quantum approach (solid 
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electronic states. 
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FIG. 12: (a) Population differences between different subband 
levels, as wells as total electronic sheet density as a function 
of the temperature, as computed from the quantum model, 
(b) Comparison between the effective plasma frequencies of 
the two resonances in the case without coupling (dashed lines) 
and when the plasmon-plasmon coupling is taken into account 
(solid lines). 



temperatures: 7 K, 30 K, 50 K and 100 K. The cavi- 
ties used in the experiments had a patch width ranging 
from 6 to 21 jam. This permits to explore a frequency 
range from 1.8 THz to 6 THz for the fundamental mode. 
The comparison between the experimental spectra (dot- 



ted lines) and the effective medium electromagnetic simu- 
lations (continuous lines) are presented in Fig (13] at T =7 
K. For simplicity, we consider only measurements with an 
incident angle of ^ = 10°, where only the odd resonances 
are excited. The agreement between the experiments and 
electromagnetic simulations is excellent. The simulations 
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FIG. 13: (a) and (b) Experimental reflectivity spectra at 
= 10° (dotted curves) and simulations (solid lines) for all 
measured cavities, at low temperature T —7 K. The number 
next to each spectra indicates the size of the patch s. 



reproduce correctly both the positions of the different 
resonances and the evolution of their reflectivity dips. 

To summarize the comparison between measurements 
and experiments in Fig(ll] we have plotted the polari- 
ton dispersion curves for the four different temperatures 
above mentioned. Again, very good agreement is ob- 
tained between the experiment (triangles) and the effec- 
tive medium model based upon the dielectric constant 
Eq. (p!2]) (continuous lines) for every temperature. These 
dispersion curves indicate clearly how the splitting of 
the high energy excitation 2^3 increases while that 
of the 1^2 excitation decreases as the temperature 
is raised. In accordance the polariton gap associated 
with the 1^2 excitations progressively decreases with 
the temperature and at T =100 K no gap can be seen. 
The "temperature knob" therefore allows us to tune the 
system from the ultra-strong coupling to the "ordinary" 
strong-couphng regime [7] . 

The study of the polaritonic dispersion brings insights 
on the impact of the dissipation (linewidth enlargement) 
on the measured splitting. In a system with loss the 
observable Rabi splitting is in general smaller than the 
actual strength of the light-matter interaction |13i]. For 
an ideal lossless system, the polariton dispersion can be 
obtained analytically from the eigenvalue equation of the 
quantum Hamiltonian [10|]. For our particular system, 
the dispersion equation for two plasmon modes coupled 
with a lossless TMq mode is: 



1 



P2 



P2 



wi 



(13) 



The real roots of Eq. (p!3|) are plotted in a dashed-line 
in FigdSl for T =50 K, where both ISB excitations are 
clearly present. The results from the dielectric constant 
model, that fits exactly the measured data, are indicated 
as a solid-line. Not surprisingly the predicted Rabi split- 
tings for the ideal system described by Eq. (p!3|) are greater 
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FIG. 14: Dispersion curves for the coupling between the two 
ISB plasmons with the fundamental cavity mode. The solid 
lines are electromagnetic simulations, and the triangles are 
measured polariton frequencies, (a) T =7 K, (b) T =30 K, 
(c) T =50 K, (d) T =100 K. The grey area indicates the 
polariton gap for the 1^2 ISB plasmon. 



than those obtained experimentally. This is the results 
of the reduced coherence in the real system due to the 
loss introduced by the metallic layers of the cavity, the 
finite width of the ISB resonance, as well as the decreased 
transparency of the semiconductor close to the restrahlen 
band. 

For an infinitely coherent system, according to Eq. (p!3|) 
the polariton splittings 20.^ and 21^23^ would be ex- 
actly Ri^2 and i^2^3- Interestingly, a very good fit will 
be obtained also by replacing in Eq. (p!3|) Ri^2 and i^2^3 
with the experimental values for the splittings (2Qf^ and 
21^23^) shown in the dotted-curve of FigdSl This 

implicates that the polaritonic dispersion can be very well 
reproduced from the perfect quantum model, as long as 
the measured values of the splittings 2^1^^^ are used for 
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FIG. 15: Dispersion curves computed from the different mod- 
els mentioned in the text. 



TABLE I: 1 ^ 2 transition 



T(K) 


Rl^2 " 


^ 




-L cav 


Ci^2 


7 


1.64 


1.40 


0.20 


0.55 


12.77 


30 


1.41 


1.20 


0.53 


0.55 


3.4 


50 


1.23 


1.04 


0.60 


0.55 


2.3 


100 


0.90 


0.70 


0.65 


0.55 


1.1 



^In THz, from the effective medium model. 
^In THz, from the dispersion curves. 

^In THz, simulation of a s =12.5/im cavity without carriers. 



the light-matter coupling constants in the Hamiltonian. 
Moreover, it justifies the approach of several experimen- 
tal groups that fit the measured dispersion curves with 
fully coherent equations, like Eq.([T3]), by using the ex- 
perimentally derived Rabi frequency at resonance [5-8]. 

The comparison between the ideal system and our real 
system is summarized in Table |T] and Table |TT1 for the 
1^2 transition and 2^3 transition, respectively. A 
sizeable discrepancy can be measured between the mod- 
elled coupling constant i^2^3 and the observable split- 
ting, 2f]23^, particularly when the linewidth is an impor- 
tant fraction of these quantities. To quantify this effect 
we define cooper at ivity parameter, as in atomic physiscs, 
that sets the onset for the strong coupling in our sys- 
tem. In atomic physics, where an atom is interacting 
with a microcavity mode, the cooperativity parameter is 



TABLE II: 2 ^ 3 transition 



T(K) 




ryr\exp h 




-L cav 


^^2^3 


7 


0.72 


0.24 


0.48 


0.77 


0.7 


30 


1.18 


0.61 


0.58 


0.77 


1.6 


50 


1.46 


0.80 


0.68 


0.77 


2.0 


100 


1.35 


0.82 


0.70 


0.77 


1.7 



^In THz, from the effective medium model. 
^In THz, from the dispersion curves. 

•^In THz, simulation of a s =7.8/im cavity without carriers. 



defined as the ratio between the rate of photon scatter- 
ing into the microcavity mode divided by the rate of the 
free-space scattering (radiative enlargement of the atomic 
transition) [26, 27]. This definition can be adapted to our 
system by taking into account the fact that the dominant 
contribution to linewidth enlargement is not of radiative 
origin. Indeed, the decay rate of the ISB plasmon is rj/2 
and is determined mainly by non-radiative mechanisms, 
such as phonon relaxation and interface scattering. We 
can define thus the cooperativity parameter Cj as: 



where R^j/AV^q^^ is the rate of photon emission in the 
cavity mode as provided by Fermi Golden rule [28], since 
the coupling constant of the Hamiltonian of Eq. (p!Q|) is 
Rj /2 and Fcav is the full- width at half maximum of the 
cavity mode. With this definition, the criterion for the 
intersubband system to be in the strong coupling regime 
is Cj > 1. 

The values of Cj are reported in the rightmost column 
of the Tables |I] and [III and show that this definition is 
in accordance with our data. In particular we can ob- 
serve that the coupling between the 2^3 transition and 
the fundamental mode has in general a lower coopera- 
tivity and therefore a much higher discrepancy between 
the values of the coupling constant Rj and the measured 
splittings 2(^23^ is observed. This is due to the dissipa- 
tions effect of the 2^3 transition being more impor- 
tant than those of the 1^2 transition, probably due 
to the proximity of the GaAs phonon absorption band. 
According to the criterion Cj > 1, the 2^3 transi- 
tion is not in strong coupling at T =7 K as it can be 
also seen in Fig JIOT c) and Fig lHT a). More generally, the 
importance of the dissipation effects on the onset of the 
strong coupling could explain the discrepancies between 
the light-matter coupling estimated from simple analyt- 
ical formulas like Eq.([T]) and experimental observations, 
especially in a case of large ISB linewidths or of metallic 
cavities with high losses as those reported in Ref.^8|]. 

IV. CONCLUSION 

In conclusion, we have presented a study of intersub- 
band systems strongly coupled with CD micro-resonators. 
Our analysis permits to relate the polariton spectroscopic 
features to microscopic details of the heterostructure, 
such as the thermal distribution of electrons on the dif- 
ferent subbands and the coupling between different col- 
lective oscillations in the quantum well. We interpret 
our experimental results using an effective medium ap- 
proach, which connects the quantum simulation of the 
heterostructure with the electromagnetic modelhng of 
the microresontor containing the quantum well slab. We 
also discuss some of the important aspects of the experi- 
mental system, such as finite linewidths or charge transfer 
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occurring at metal- semiconductor interfaces. Such mod- 
elling can be further refined by including, for instance, 
inhomogeneous broadening effects owed to disorder like 
in the numerical approach proposed in Ref.[25]. Since we 
consider samples with high electronic densities, we expect 
that other fundamental effects, like exchange-correlation 
and excitonic like shift [29] to be negligible. 

We therefore believe that the approach described 
here, which is a combination of microscopic self- 
consistent modelling of the heterostructure and electro- 
magnetic simulations through the effective medium con- 
stant fEq. (p!2]) ) could be applied as a designer tool for 
ISB light-matter coupled quantum devices. 
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Appendix A: Periodic Poisson-Schrodinger solver 

In a typical self-consistent Poisson-Schrodinger prob- 
lem one determines the Hartree correction of the poten- 
tial Vh{z)^ that includes the contribution of the static 
charges through the Poisson equation [2|]. This contri- 
bution is added to the heterostructure potential, and 
the Schrodinger equation is solved numerically, typically 
trough the shooting method. 

As seen from FiglHl the doped quantum wells cre- 
ate a periodic potential. Hence the Hartee correction 
Vh{z) to the heterostructure potential, that depends on 
the wavefunctions is also periodic. However, such self- 
consistent Poisson-Schrodinger problem is inconvenient 
to treat with the shooting method, since the energy levels 
of the quantum wells are heavily degenerated. We there- 
fore treat only a single period with extension d (such as 
Vh{z d) = Vh{z)). The most general solution for the 
Poisson equation for a 2D charge density p{z) is then 



Vniz) = — 
eeo L 



z')p{z')dz' + ^ 



z'p{z')dz' 



(A3) 

It is easy to verify with this formula that both the 
potential Vh{z) and its derivative dVH{z)/dz are d- 
periodic. Hence Eq. (|A3[) describes the unique solution 
of the Poisson equation on a globally neutral cylindrical 
space with a circumference d. 



Appendix B: Multiband Hamiltonian 

The diagonalization of the Hamiltonian ([7|) is per- 
formed in two steps. First, we diagonalize the matter 
part of the Hamiltonian which includes the inter-plasmon 
interactions. The result is a set of new independent nor- 
mal collective plasmon modes, which are coherent super- 
positions of plasmons from different transitions. We then 
explicit the interaction with the normal modes with the 
light mode. The new coupling constants of the light- 
matter interaction then describe the collective oscillator 
strength of the normal modes. 

We therefore start with the matter Hamiltonian: 



a 

+ ^ l^^a^ipi -^Pa){pl^P^) (Bl) 

If there are N optically active plasmons in the system, 
then, after the diagolization is performed, there will be 
N independent normal modes. The bosonic operators Pj 
of the normal modes are linear combination of the initial 



ones: 



Pj = Y.{YjaPa^Tjapl) 



(B2) 



Vh{z) = — [\z- z')p{z')dz' ^Az^B 
££o Jo 



(Al) 



Here p{z) contains both the contribution from the 
wavefunctions and the donor impurities. We choose 
B = and we require that the charge density is also 
periodic: p{z -\- d) = p{z). Moreover, we require that the 
starting point z = of the period is chosen so that the 
period is globally neutral: 



Jo 



p{z')dz' = 



(A2) 



Hence the constant A is uniquely determined and 
Eq. ()Aip becomes: 



The latter are eigenmodes of the Hamiltonian ()B1 
with eigenfrequencies Wj only if: 



hWjPj 



(B3) 



Using the symmetry property ^^(3 = 'BfBa and the re- 
lation: 



[Pa^Hmat] = ^aPa + ^ K^a^{p\ + P^) (B4) 



as well as its hermitian conjugate, Eq.(jB3]) leads to the 
following eigenvalue problem: 



MVj = Wj\j 



(B5) 
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Here M is the 2N x 2N Hopfield-Bogoliubov matrix 
of the problem, and and Vj is the 2N rank eigenvector 
which corresponds to the eigenvalue Wj: 



M 



/ II Ci2 
Ci2 I2 



UJoc 



ClAT 
^2N 



\ CiAT C2Ar • • • Iat / 



(B6) 



(B7) 



Therefore the equation providing the light-coupled po- 
lariton states is: 



CO 



E 



P2 , ,2 



(B13) 



In the case where there are Nqyj initial plasmons 
are uncoupled, but oscillate at identical frequencies, the 
above formalism predicts that there is one bright excita- 
tion, with a collective oscillator strength that is ^J Nqyj 
higher. This is easily seen from the above equation, by 
setting Wj = const. 

The effective dielectric constant of the system is then 
defined through the relation eefi{oj)oj'^ = suol^^: 



Yj =^ {YjuTj^,---Yjm,Tjm) 



(B8) 



Since the matrix M real and block-symmetric, it al- 
lows reel eigenvalues with reel eigenvectors. Moreover, if 
Wj is an eigenvalue, then —Wj is also eigenvalue of M. 
Therefore the matrix M has N distinct positive eigenval- 
ues Wj which are the frequencies of the normal plasmon 
modes. The latter can be easily inferred from the numer- 
ical solution of the linear problem described by Eq. ()B5p . 

Let us now explicit the coupling of the normal modes 
with the light mode. Since the eigenvectors Vj are real, 
then it is easy to show that we have: 



Ja 



Tj^){Pa+pl) 



(B9) 



This relation, seen as a matrix- vector product, allows 
to determine + Pa as a function of Pj + Pj. Let us 
define the N x N matrix inverse: 



^aJ = {yja + Tja) 



(BIO) 



This leads to the light-matter coupling Hamiltonian of 
Eq. (p!Q|) . where the effective plasma frequencies Rj of the 
normal modes are defined as: 



eeff(^) 



E 

J 



(B14) 



The treatment of the problem is so far energy con- 
serving. When introducing Markovian relaxation mecha- 
nism with a rate rj/2 and considering a basis of coherent 
(semiclassical) states the above equation is transformed 
into Eq.([l2|) [30]. 

Let us discuss briefiy the case of a heterostructure con- 
sisiting of A^qw identical quantum wells. For simplicity, 
let us consider only two plasmons from each quantum 
well, coupled by a constant S^b, and we consider that 
the plasmons from different wells are uncoupled, owe to 
a vanishing overlap of the wavefunctions between differ- 
ent wells. In this case the matter Hamiltonian writes: 



2=1 i=l 
^ Nqw 

-^-hEab {pU^PaMvU (B15) 

The light matter-part interaction Hamiltonian is re- 
spectivelly 



Rj 



Wj 



fo fw 



(Bll) 



At this point, the Hamiltonian ([TO]) describes N inde- 
pendent quantum oscillators coupled with a single cavity 
mode. Its Hopfield matrix is very similar to Eq. ()B6p . 
except that the coupling terms are non zero only on the 
first two columns and first two rows. The characteristic 
polynomial of this problem can then be computed explic- 
itly, and the result is: 



A^QW 

Hint = ihQa Y '^^^ ~ ^Pa,i) 
Nqw 

^ihnb (^^ - ^)ipli^Pb,i) (BI6) 

i=l 

To simplifiy these Hamiltonians, we introduce the fol- 
lowing new bosonic operators [^: 



w 



of 



E 



2, ,2 



w] 



} (B12) 



Nqw 



r^= P{a,b),i 

Vqw . 1 



(B17) 
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Here the normalization factor l/y^Nqw arises from 
the requirement that [7ra,6,7r^5] = 1. Let us then com- 
pute the commutator of ttJ with the Hamiltonian ()B15p : 

[Hmat.Trl] = KJja r- ^ 

V^QW 

K:Ja7rl + ^hEabilTb + TT^^) (B18) 

Similar commutation relation is valid for tt^. These 
commutation relations are identical if we replace Eq. 
(|B15p with the effective Hamiltonian: 



^^hEab{7rl + 7ra)(7r^^ + TTb) (B19) 

Expressing the interaction Hamiltonian with the new 
operators we obtain: 

Hint = ihQa^/N^ia"^ - a)(7rj + iTa) 

^mnby^NQwia"^ - a){7rl + tt^) (B20) 

Formally, the resulting Hamiltonian is equivalent to 
Eq.©, except that the light-matter coupling constants 
for a single quantum well are multiplied by y^A^qw- This 
justifies our approach in part IHII 
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